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I.  INTRODUCTION 


Many  projectiles  used  by  the  Army  are  slender,  spin  stabilized 
bodies  of  revolution.  The  distorted  boundary  layer,  which  results  from 
a spinning  body  at  angle  of  yaw,  generates  Magnus  forces  and  moments. 
Recent  Army  interest  in  achieving  increased  range  and  greater  payload 
capacity  in  artillery  projectiles  has  led  to  designs  with  long,  slender 
ogives,  increased  projectile  length,  and  boattailed  afterbodies.  These 
designs  have  resulted  in  decreased  drag  and  an  increase  in  range; 
however,  the  gyroscopic  stability  of  these  shapes  is  less  than  that  of 
more  conventional  designs.  This  means  that  these  new  shapes  are  more 
susceptible  to  a Magnus  induced  instability.  The  Magnus  force  is 
small,  typically  1/10  to  1/100  of  the  normal  force;  however,  its  effect 
is  important  because  the  Magnus  moment  acts  to  undamp  the  projectile 
throughout  its  flight.  Thus,  it  is  desirable  to  minimize  the  Magnus 
moment  so  that  the  projectile  flies  at  a small  average  angle  of  attack 
and  achieves  the  greatest  range  capability. 

BRL  has  developed  numerical  techniques  for  computing  Magnus 
effects  (forces  and  moments)  which  result  from  spin  induced  distortion 
of  the  boundary  layer.  The  computational  procedure  is  described  in 
references  1,  2,  3,  and  4 and  consists  of  numerical  techniques  for 
computing:  (1)  the  three  dimensional  supersonic  inviscid  flow  over  a 
yawed  body  of  revolution;  (2)  the  three  dimensional  laminar/turbulent 
boundary  layer  development  over  a yawed,  spinning  body  of  revolution; 

(3)  the  three  dimensional  boundar>  layer  displacement  surface;  and 

(4)  the  three  dimensional  supersonic  flow  over  a yawed  body  plus  three 
dimensional  boundary  layer  displacement  surface  to  yield  pitch  and  yaw 


J.  H.  A.  Duyer^  "Three  Dimensional  Flau  Studies  Over  a Spinning  Cone 
at  Angle  of  Attaaky"  BRL  Contract  Report  No,  12?,  February  1974, 
U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland.  AD  774795. 

2,  H.  A.  Dwyer  and  B.  R.  Sanders,  "Magnus  Forces  on  Spinning  Super- 
sonic Cones.  Part  I:  The  Boundary  Layer, " BRL  Contract  Report  No. 
248,  July  1975,  U.S.  Army  Ballistic  Research  Laboratory , Aberdeen 
Proving  Ground,  Maryland.  AD  A012518.  Also  AIAA  Journal,  Vol.  14, 
No.  4,  April  1976,  p.  498, 

2,  B.  R,  Sanders,  "Three-Dimensional,  Steady,  Inviscid  Flaw  Field 

Calculations  With  Application  to  the  Magnus  Problem,"  PhD  Disserta- 
tion, University  of  California,  Davis,  California,  May  1974. 


4,  W.  B.  Sturek  (et  al),  "Computations  of  Turbulent  Boundary  Layer 
Development  Over  a Yawed,  Spinning  Body  of  Revolution  With  Appli- 
cation to  Magnus  Effect,"  BRL  Report  No.  1985,  May  1977,  U.S.  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland. 
AD  A041228. 
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plane  aerodynamic  coefficients,  Experimental  studies  are  being  carried 
out  to  provide  data  which  will  help  guide  this  theoretical  effort.  The 
purpose  of  the  experimental  measurements  reported  here  is  to  provide 
data  to  evaluate  the  theoretical  computations  of  turbulent  boundary 
layer  development.  The  objective  of  this  report  is  to  provide  a summary 
of  the  boundary  layer  profile  data  with  some  analysis  and  comparison 
with  theory.  A complete  tabulation  of  experimental  data  is  presented  in 
Part  II  of  this  report  (reference  5). 


II.  EXPERIMENT 

All  experimental  results  presented  in  this  report  were  obtained 
on  the  secant-ogive-cylinder  model  (SOC)  shown  in  Figure  1.  The  model 
is  57.15  mm  (2.25  inches)  in  diameter  and  342.9mm  (13.5  inches)  long. 

A boundary  layer  trip  was  placed  on  the  ogive  to  insure  the  location  of 
the  start  of  turbulent  flow.  The  tests  were  conducted  in  the  BRL 
Supersonic  Wind  Tunnel  No.  1 which  is  a continuous  flow  tunnel  with  a 
test  section  of  330  x 381  mm  (13  x IS  inches).  Measurements  of  the 
total  head  pressure  through  the  boundary  layer  were  made  with  a flattened 
impact  pressure  probe  1.5  mm  wide  by  0.15  mm  high.  The  probe  was 
electrically  isolated  from  the  probe  holder  so  that  contact  with  the 
model,  for  non-spinning  runs,  could  be  determined  with  an  ohmmeter. 

The  probe  drive  mechanism  moved  the  probe  perpendicular  to  the  model 
centerline.  In  addition,  the  probe  drive  mechanism  could  be  positioned 
circumferentially  about  the  model.  Figure  2 shows  the  model,  probe,  and 
probe  drive  unit  installed  in  the  tunnel.  Figure  3 is  a shadowgraph 
showing  the  flow  over  the  model  and  probe  for  4.2  degrees  angle  of 
attack. 

The  boundary  layer  survey  procedure  was  to  bring  the  impact  probe 
from  outside  the  boundary  layer  down  toward  the  model  and  touching  it 
for  the  no-spin  case.  Immediately  following  a no- spin  run,  the  model 
was  brought  up  to  the  333  rps  (20,000  rpm)  spin  rate  and  the  probe  was 
brought  down  through  the  boundary  layer  to  within  approximately  0.1  mm 
from  the  surface.  Data  were  obtained  at  3.33,  4.44,  and  5.56  calibers 
(body  diameters)  from  the  nose  and  at  angles  of  attack  from  0 to  6.3 
degrees.  Data  were  acquired  circumferentially  in  30  degree  increments 
and  also  at  10  degrees  on  each  side  of  the  leeward  ray  (4i  = 180  degrees). 
At  most  positions,  surveys  were  made  at  both  0 and  333  rps:  the  spin 
rate  of  333  ips  corresponds  to  a dimensionless  spin  rate  (pd/V)  of  0.19 
at  Mach  3.0.  Tunnel  conditions  for  the  tests  were:  Mach  3.0;  a supply 
temperature  of  310  K;  a supply  pressure  of  298  kPa.  These  conditions 
provided  a Reynolds  number  of  7.3  x 10^  based  on  model  length.  Local 
Mach  numbers  within  the  boundary  layer  were  determined  from  the 


5.  L,  D.  Kayser  and  W,  B.  Sturek,  "Experimental  Measurements  in  the 
Turbulent  Boundary  Layer  of  a Yauied,  Spinning  Ogive-Cylinder  Body 
of  Revolution  at  Maah  3.0.  Part  II.  Data  Tabulation,"  to  be 
published  as  a BEL  Memorandum  Report. 
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Rayleigh  pitot  formula  assuming  a constant  static  pressure  across  the 
boundary  layer.  The  data  in  this  report  were  reduced  using  the  experi- 
mental values  of  wall  static  pressure  obtained  by  Reklis®.  Figure  4 
is  a comparison  of  experimental  surface  pressures  and  theoretical 
surface  pressures  computed  with  the  inviscid  program  discussed  in  the 
introduction.  The  model  surface  temperature  was  assumed  to  be  equal  to 
the  adiabatic  wall  temperature  for  turbulent  flow  and  the  recovery 
factor  was  taken  as  the  cube  root  of  the  Prandtl  number.  The  temperature 
distribution  in  the  boundary  layer  was  found  by  assuming  the  Crocco 
linear  total  temperature-velocity  relationship; 


With  temperature,  pressure,  and  Mach  number  known,  local  values  for 
density  and  velocity  have  been  computed  and  tabulated.  The  integral 
parameters — displacement  thickness,  momentum  thickness,  and  velocity 
thickness — have  been  computed  by  integrating  the  profile  data. 

The  probe  axis  was  aligned  longitudinally  with  the  model  axis. 

Some  uncertainty  is  inherent  in  the  profile  data  due  to  the  probe  not 
being  aligned  with  the  local  flow  direction  within  the  boundary  layer. 
The  uncertainty  due  to  cross  flow  would  be  of  the  order  of  angle  of 
attack  at  the  outer  edge  of  the  boundary  layer  when  probing  the  sides 
of  the  model  (4)  = 90  and  270  degrees) . The  uncertainty  due  to  the 
effects  of  spin  on  the  crossflow  velocity  are  expected  to  be  small 
because  the  large  velocity  gradients  in  a turbulent  boundary  layer  would 
confine  the  greatest  effect  of  flow  angularity  to  a very  small  region 
near  the  surface  which  cannot  be  probed  accurately  using  a total  head 
probe. 

Wall  shear  stress  was  obtained  for  the  non-spinning  model  using 
the  Preston  tube  technique.  The  Preston  tube  is  a circular  total  head 
probe  which  is  brought  down  to  the  model  surface  and  is  designed  in 
size  to  lie  within  the  logarithmic  portion  of  the  law-of-the-wall 
velocity  profile.  The  wall  shear  stresses  are  then  computed  from  the 
Preston  tube  measurements,  the  Preston  tube  size,  and  model  surface 
pressures  using  the  correlation  relations  found  in  reference  7.  The 
free  stream  Mach  number  of  3.0  was  used  in  reducing  the  Preston  tube 
data  rather  than  the  local  Mach  number  at  the  survey  station. 

6.  P.  P.  Reklia  and  W,  B,  Sturek^  "Meaaurements  of  Wall  Static 
Preaaure  on  Slender  Bodies  of  Revolution  at  Angle  of  Attaakf" 
to  be  publiahed  as  a BRL  Memorandum  Report. 

7.  P.  Bradshau  and  K.  Unsworth^  "A  Rote  on  Preston  Tube  Calibrations 
in  Compressible  Flow,"  IC  Aero  Report  72-07,  September  1972, 

Imperial  College  of  Science  and  Technology , London,  Great  Britain. 


I 


9 


Data  Accuracy  — The  pressure  transducers  are  linear  to  within 
± 0.25%  of  full  scale  value.  The  data  acquisition  system  measurement 
accuracy  is  approximately  0.1%  of  full  scale.  Full  scale  is  rarely 
achieved  on  the  transducer  or  the  system  measurement  range;  therefore, 
the  accuracy  of  measured  pressures  is  estimated  to  be  ± 1.0%. 
Contributing  errors  in  determining  the  probe  height  relative  to  the 
model  surface  are  as  follows:  (1)  accuracy  of  probe  height  calibration; 
(2)  measurement  error;  and  (3)  error  in  determining  exactly  when  the 
probe  contacted  the  model.  The  overall  accuracy  of  the  probe  height 
value  (y-coordinate)  is  estimated  to  be  within  ±0,1  mm.  Values  of 
skin  friction  obtained  by  the  Preston  tube  technique  are  estimated  to 
be  accurate  to  within  ± 15%. 


III.  DISCUSSION  OF  RESULTS 

A summary  of  test  conditions  for  the  SOC  boundary  layer  survey  data 
is  given  in  Table  I.  A complete  set  of  boundary  layer  profile  data  and 
Preston  tube  skin  friction  data  are  presented  in  Part  II  (reference  5) 
of  this  report. 

To  help  clarify  the  data,  the  orientation  of  the  probe  with 
respect  to  the  model  must  be  known.  The  circumferential,  or  azimuthal, 
position  on  the  model  of  4)  = 0 degrees  is  defined  as  the  most  windward 
ray  on  the  model  when  the  model  is  at  some  angle  of  attack.  Looking 
upstream  at  the  model  base,  with  the  model  at  positive  angle  of  attack, 

4 = 0 is  on  the  bottom  (6  o'clock);  (fi  = 90  degrees  is  to  the  left 
(9  o'clock);  (j)  = 180  degrees  is  on  the  top  (12  o'clock);  and  4 = 270 
degrees  is  to  the  right  (3  o'clock).  A clockwise  spin  is  positive; 
therefore,  a positive  spin  gives  a surface  velocity  in  the  same 
direction  as  cross  flow  on  the  left  side  of  the  model.  On  the  right 
side,  cross  flow  and  model  surface  velocities  are  in  the  opposite 
direction.  The  model  coordinate  system  is  shown  in  Figure  5 with  the 
arrows  indicating  positive  directions. 

Figure  6 is  a comparison  of  theoretical  and  experimental  velocity 
profiles  at  zero  spin  and  4.2  degrees  angle  of  attack.  This  figure 
illustrates  the  thickening  of  the  boundary  layer  when  moving  from  the 
windward  side,  4 = 0,  to  the  leeward  side,  (p  = 180.  On  the  leeward 
side,  the  theoretical  velocities  are  greater  than  experimental  values 
near  the  model  surface  eind  the  theoretical  velocities  are  smaller  than 
experimental  values  near  the  edge  of  the  boundary  layer.  On  the  wind- 
ward side,  though  not  very  noticeable,  the  situation  is  reversed.  The 
differences  in  profile  shape  will  give  compensating  effects  when  compu- 
ting integral  parameters.  Figure  7 is  a comparison  of  theoretical  and 
experimental  velocity  profiles  for  the  spinning  model  case.  The 
differences  between  theory  and  experiment  are  virtually  the  same  as 
for  the  no-spin  case  of  Figure  6.  The  effect  of  spin  on  experimental 
velocity  profiles  is  shown  in  Figure  8 where  profiles  on  the  left  side 
of  the  model  are  compared  with  those  on  the  right  side.  On  the  wind- 
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ward  side  of  the  model  (ij!  = 0 to  90  and  360  to  270),  there  is  almost 
no  measurable  effect  of  spin.  However,  on  the  leeward  side  at  <ti  = 120 
vs  240  and  ip  = 150  vs  210,  the  profile  shapes  differ  substantially.  The 
effect  of  cross  flow  in  opposition  (cj)  = 180  to  360)  to  surface  model 
rotation  (e.g.  ip  = 210)  is  to  decrease  the  fullness  of  the  profile 
which,  of  course,  will  result  in  a larger  displacement  thickness.  It 
is  also  of  interest  to  note  that  the  primary  effect  of  spin  is  to 
change  the  profile  shape  rather  than  to  change  the  total  thickness. 

Values  for  tiie  longitudinal  component  of  displacement  thickness 
are  compared  in  Figure  9.  The  agreement  between  theory  and  experiment 
is  generally  good;  however,  it  is  seen  at  the  forward  station,  4>  = 0, 
(2/D  = 3.33)  that  theoretical  thicknesses  are  slightly  greater  than 
experiment  and  at  the  aft  station  theoretical  thicknesses  are  slightly 
smaller  than  experiment.  This  situation  indicates  that  the  boundary 
layer  actually  grows  at  a faster  rate  than  predicted  by  theory;  however, 
this  is  not  particularly  surprising  since  the  turbulence  model  did  not 
provide  for  any  adjustment  as  a function  of  pressure  gradient.  The 
effect  of  spin  on  displacement  thickness  6*^  can  be  seen  in  Figure  10 

where  the  increment  of  6*^  due  to  spin  is  plotted  on  an  expanded  scale 

for  a = 4.2  degrees.  The  effect  on  displacement  thickness  is  seen  to 
be  significant  only  in  the  vicinity  of  the  leeward  side  ((f  = 180 
degrees) . The  agreement  between  theory  and  .experiment  is  encouraging 
evidence  that  the  numerical  technique  accurately  models  the  effect  of 
surface  spin. 

Measured  values  for  skin  friction  coefficient  obtained  using  the 
Preston  tube  technique  are  compared  to  theory  in  Figure  11.  The  skin 
friction  coefficient  is  referenced  to  free-stream  static  properties 
upstream  of  the  model  rather  than  the  more  conventional  approach  of 
using  local  properties  at  the  edge  of  the  boundary  layer.  The  agreement 
indicated  is  within  ± 10%.  This  is  considered  quite  good  since  the 
Preston  tube  is  expected  to  yield  an  accuracy  of  ± 10%  for  two  dimen- 
sional flat  plate  boundary  layer  flow.  The  use  of  the  Preston  tube  to 
obtain  measurements  in  a three  dimensional  boundary  layer  flow  using 
two  dimensional  calibration  data  must  be  regarded  as  speculative  and 
mainly  of  qualitative  interest. 

Experimental  displacement  thickness  data  at  6.34  degrees  angle  of 
attack  are  shown  in  Figure  12.  The  significant  difference  from  the  4.2- 
degree  case  is  the  dip,  or  decrease,  in  6*^  near  the  leeward  ray  ($  = 

180  degrees).  The  thicker  boundary  layers  (15*^^)  on  either  side  of  the 

leeward  ray  are  believed  to  be  caused  by  the  existence  of  longitudinal 
separation  type  vortices  which  are  beginning  to  develop.  Such  vortices 
could  create  local  areas  of  favorable  and  adverse  pressure  gradients 
that  would  cause  the  complexity  illustrated  in  Figure  12.  This 
phenomenon  is  evidence  that  conventional  boundary  layer  theory  will  not 
be  adequate  at  the  higher  angles  of  attack. 
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IV.  CONCLUDING  REMARKS 
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The  primary  cc/ntribution  of  the  experiment  described  in  this 
report  is  the  unique  set  of  data  which  can  be  used  to  evaluate 
theoretical  computational  procedures  for  three-dimensional  turbulent 
boundary  layer  development.  Analysis  of  these  data  and  comparisons  to 
theory  support  the  following  statements  which  summarize  the  primary 
findings  of  this  report. 

(1)  The  boundary  layer  data  show  the  effect  of  spin  to  be 
confined  primarily  to  the  leeward  side  of  the  model. 

(2)  The  excellent  agreement  between  theory  and  experiment  for 

is  encouraging  evidence  that  the  boundary  layer  theory  accurately 

models  the  effects  of  surface  spin. 

(3)  The  circumferential  and  longitudinal  growth  of  the  boundary 
layer  show  reasonably  good  agreement  between  theory  and  experiment. 
However,  there  is  disagreement  in  the  longitudinal  growth  rate  of  the 
boundary  layer  between  theory  and  experiment. 

(4)  There  is  a consistent  tendency  for  the  theory  to  predict  a 
velocity  profile  that  is  more  full  than  experiment  near  the  wall  and 
less  full  than  experiment  near  the  edge  of  the  boundary  layer. 

(5)  If  the  dip  in  the  curve  of  Figure  12  is  an  indication  of 
longitudinal  development  of  vortices,  significant  flow  components  in 
the  y-direction  would  exist.  Conventional  boundary  layer  equations 
would  not  be  expected  to  be  valid  beyond  approximately  5-degrees  angle 
of  attack  since  the  equations  do  not  include  y-momentum. 
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Figure  2.  Model  Installation  Photograph 


Figure  3.  Shadowgraph,  a = 4.2  degrees,  Z/D  = 3.33 


Figure  6.  Velocity  Profiles,  Theory  Compared  With  Experiment 


files.  Effects  of  Spin 


Table  1.  Test  Run  Summary 


Z/D  = 3.33 


i 

a = 0 

a = 2.1 

a = 4.2 

a = 5.3 

a = 6,. 

0 

0 

0,333* 

0,333 

0,333 

30 

0 

0,333 

0,333 

333 

60 

0 

0,333 

0,333 

0,333 

90 

0 

0,333 

0,333 

120 

0,333 

0,333 

0,333 

0,333 

150 

0,333 

0,333 

0,333 

0,333 

170 

0,333 

0,333 

0,333 

0,333 

180 

0,333 

0,333 

0,333 

0,333 

190 

0,333 

0,333 

0,333 

0,333 

210 

0,333 

0,333 

0,333 

0,333 

240 

0,333 

0,333 

0,333 

0,333 

270 

0,333 

0,333 

300 

0,333 

0,333 

0,333 

330 

0,333 

0,333 

0,333 

Z/D  = 4.44 

0 

0 

0,333 

0,333 

0,333 

30 

0 

0,333 

0,333 

0,333 

60 

0 

0,333 

0,333 

0,333 

90 

0 

0,333 

0,333 

120 

0,333 

0,333 

0,333 

0,333 

150 

0,333 

0,333 

0,333 

0,333 

170 

0,333 

0,333 

0,333 

0,333 

180 

0,333 

0,333 

0,333 

0,333 

190 

0,333 

0,333 

0,333 

0,333 

210 

0 

0,333 

0,333 

0,333 

240 

0,333 

0,333 

0,333 

0,333 

270 

0,333 

0,333 

300 

0,333 

0,333 

0,333 

330 

0,333 

0,333 

0 
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Table  1.  Test  Run  Summary  (Continued) 


i 

a = 0 

a = 2.1 

a = 4.2 

a = 5.3 

a = 6. 

Z/D  = 5.56 

0 

0 

0,333* 

0,333 

0,333 

30 

0 

0,333 

0,333 

0,333 

60 

0 

0,333 

0,333 

0,333 

90 

0 

0,333 

0,333 

120 

0,333 

0,333 

0,333 

150 

0,333 

0,333 

0,333 

0,333 

170 

0,333 

0,333 

0,333 

0,333 

180 

0,333 

0,333 

0,333 

190 

0,333 

0,333 

0,333 

0,333 

210 

0,333 

0,333 

0,333 

0,333 

240 

0,333 

0,333 

0,333 

0,333 

270 

0 

0,333 

300 

0,333 

0,333 

0,333 

330 

0,333 

0,333 

0,333 

* Model  Spin  Rate,  0 and  333  rps,  333  rps  = 20,000  rpm 
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LIST  OF  SYMBOLS 

skin  friction  coefficient,  x /q 

diameter  of  model  base,  mm 
model  spin  rate,  radians/second 
impact  probe  pressure,  kPa 

model  wall  static  pressure,  kPa 

free-stream  static  pressure,  kPa 

free-stream  dynamic  pressure,  kPa 

secant-ogive-cylinder 
tunnel  total  temperature,  K 

local  total  temperature,  K 

model  wall  temperature,  K 

velocities  in  boundary  layer  coordinates,  m/s 
velocity  at  edge  of  boundary  layer,  m/s 

velocity  along  model  trajectory,  m/s 
boundary  layer  coordinates.  Figure  5 
longitudinal  model  axis  coordinate,  mm 
angle  of  attack,  degrees 

longitudinal  component  of  displacement  thickness,  cm 

increment  of  displacement  thickness  due  to  spin,  cm 
model  wall  shear  stress,  N/m^ 

circumferential  boundary  layer  coordinate,  degrees 
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